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From first principles calculations we investigate the electronic structure and the magnetic prop- 
erties of EuO under hydrostatic and epitaxial forces. There is a complex interdependence of the O 
2p and Eu 4/ and 5d bands on the magnetism in EuO, and decreasing lattice parameters is an ideal 
method to increase the Curie temperature, T c . Compared to hydrostatic pressure, the out-of-plane 
compensation that is available to epitaxial films influences this increase in T c , although it is mini- 
mized by the small value of poisson's ratio for EuO. We find the semiconducting gap closes at a 6% 
in-plane lattice compression for epitaxy, at which point a significant conceptual change must occur 
in the active exchange mechanisms. 



The desire to connect the electronic and magnetic 
properties of materials to enable new device character- 
istics has renewed focus on the many unique proper- 
ties of EuO, one of the initially discovered ferromagnetic 
semiconductors The divalent Eu ions in EuO possess 
a very large local moment from the half filled 4/ band 
producing a saturation magnetization of 7/ib El 12 while 
a gap of 1.2eV separates the half filled 4/ band from 
the 5d6s conduction band. [3] Low levels of electron dop- 
ing are readily achieved with oxygen vacancies, EuOi-j;, 
leading to an insulator-to-metal transition (IMT) in con- 
junction with the ferromagnetic ordering|4J and 100% 
spin polarization of the conduction electrons. [5] EuO may 
therefore be a good candidate for a spin injection mate- 
rial. Cation and anion doping has been shown to increase 
the Curie temperature (T c ) up to 170K,[51|71|S] from that 
of 69K for stoicheometric EuO .[2 

The integration of EuOi-z with Si and GaN has been 
successfully demonstrated [HE] so there is renewed inter- 
est in increasing T c towards room temperature to enable 
spintronic applications. Beyond the effects of doping on 
T c , it has been shown that a hydrostatic pressure of about 
100 kbar can increase T c to 200 K.[10] Although hydro- 
static pressure is not an option for device applications, a 
similar effect on T c may be obtainable by using epitaxial 
strain. However, epitaxy generates a biaxial stress state 
which is quite different from the isotropic stress state 
of hydrostatic pressure. Furthermore, there are a num- 
ber of competing exchange mechanisms that determine 
the magnetism in EuO, and is is not known exactly how 
this competition plays out as the lattice parameters are 
changed. 

In this paper we explore the effects of biaxial stress 
in comparison with isotropic stress on the bandstructure 
and magnetism of EuO with density functional theory 
(DFT). We find that the biaxial stress state expected in 
an epitaxial film will lead to a 50% increase in the T c 
before an insulator to metal transition occurs. This in- 
crease, although not as great as that found for isotropic 
stress, is much larger than expected for an 3-dimensional 
material under biaxial stress due to the small value of 



Poisson's ratio for EuO. From a band structure perspec- 
tive of the magnetism in EuO, we find that reducing 
lattice parameters is an ideal method of increasing T c . 
Although the closing of the semi-conducting gap, which 
we find at 5% compression for isotropic stress and 6% 
compression of biaxial stress, must lead to a conceptu- 
ally different exchange mechanisms, no sudden changes 
are seen in the magnetic properties. The combined ef- 
fects of doping and epitaxy are explored and are found 
to generate a significant increase in the mean field T c of 
up to 175K. 

The EuO band structure is calculated with the full- 
potential linearized augmented plane-wave DFT code 
WIEN2K.[11 The exchange and correlation effects are 
treated within the generalized gradient approximation 
(GGA), after Perdew et at [12] The LSDA+U method[T3] 
was used to account for strong correlations between the 
electrons in the Eu 4/ shell. Note that standard GGA or 
local spin density approximation (LSDA) predict EuO 
to be a metal whereas measurements clearly show the 
existence of conductivity gap. [3] The exchange param- 
eter (J#=0.77 eV) was chosen according to Ref. HU 
and the on-site Coulomb repulsion for the Eu 4/ or- 
bital was set to C//=8.3eV, while that of the O p orbital 
were U p = 4.6eV and 3 H =1.2 eV.[TSl[TBl[T7] With these 
values the LSDA+U band structure with the {001} an- 
tiferromagnetic spin configuration (AFMj) or the NiO- 
type {111} AFM spin configuration (AFM//) shows a 
gap of 1.2 and 1.3eV, respectively. The ferromagnetic 
(FM) spin arrangement exhibits a gap of 0.7eV (see Fig- 
ure 0. Both are consistent with the experimental optical 
abosorption gaps of 0.9 and 1.2 eV observed below and 
above the magnetic transition temperature. [HI US] Our 
calculations show a minimum in the total energy of the 
bulk EuO at a=5.1578A, in good agreement with the 
measured latticee parameter of 5.1439A.[2"U] 

The calculations were done for the isotropic stress case 
of hydrostatic pressure, where the cubic unit cell is main- 
tained, and the biaxial stress case of epitaxial film growth 
where the unit cell becomes tetragonal. For all the biaxial 
stress case calculations, the in-plane lattice parameter is 
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FIG. 1: (a) The spin resolved density of states of bulk EuO in 
the FM configuration with optimized bulk lattice parameters, 
a=5.1578A. The bands are labeled and the zero of energy is 
at the Fermi energy, Ef. 



held fixed and the minimum in the total energy is found 
as a function of the c lattice parameter. 

Since the Eu 2+ ion is in the / high spin configura- 
tion, the magnetic properties of EuO are often described 
by a Heisenberg model. According to the description of 
Mauger and Godart[3], the nearest neighbor exchange, 
Ji, is based on a 4/ electron having a virtual excitation 
to the 5c? band where it experiences an exchange interac- 
tion with the 4/ spin on a nearest neighbor. This leads to 
ferromagnetic coupling. The nature of the next nearest 
neighbor exchange constant, J2, is thought to be rather 
complex and involves various competing exchange paths 
between the O 2p, Eu 5d and 4/ orbitals. A model cal- 
culation of J2 in EuO shows that it is still ferromagnetic, 
but only about 30% the size of Ji.[3] 

In the isotropic stress case, we can extract Ji and 
J2 coupling constants from our calculations by relating 
the total energies of the various magnetic configurations 
(FM, AFM/, AFM//)[2T] from the LSDA calculations 
with the Heisenberg model. We assume that the third 
nearest neighbor coupling, J3, is very small 22, 23 . From 
these values, shown as a function of lattice parameter 
in Figure |2^, we can calculate a mean field value for 
T c = (2/3)5(5 + l)(12Ji + 6J 2 ). Our exchange parame- 
ters for the optimized bulk lattice parameter ( J\ — 0.66, 
J 2 = 0.19) are in good agreement with neutron measure- 
ments, and previous calculations. [23] [24] 

Previous experimental studies have looked at the ef- 
fect of hydrostatic pressure on T C [1(3J US] and the unit 
cell volume [26 in EuO. In Figure [2J3 we combine that 
data to plot the experimental change in T c as a function 
of lattice parameter. Included in this figure is our cal- 
culation of the mean field T c of the isotropic case as the 
lattice parameter is changed. There is quite good agree- 
ment between the experimental and calculated mean field 
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FIG. 2: (a) The nearest neighbor and next nearest neigh- 
bor exchange constants, Ji and J2, as a function of lattice 
parameter for the isotropic stress case, (b) The calculated 
mean field T c as a function of in-plane lattice parameter for 
the isotropic stress case and the biaxial stress state of an epi- 
taxial film. Both solid lines change to dashed lines at their 
respective IMT. Included in the graph is the experimental T c 
- or Ti which tracks T C [TD] - as a function of lattice parame- 
ters compiled from Ref 26 25 (circles) and 10 (triangles). The 
inset shows the change in the out-of-plane lattice parameters 
in the biaxial stress case. 



T c , with the expected over estimation of the mean field 
value. 

The calculations for the isotropic stress case indicate 
the closing of the Eu 4/-Eu 5d6s gap, or IMT, at a/ao — 
0.95. This compares well with the appearance of a drude- 
like peak in the optical reflectivity and a jump in the 
room temperature resistivity at a pressure of w 130 kbar, 
or a/a — 0.96.(101 US] As seen in Figure 2, a saturation 
of the experimental parameter Tj, which is thought to 
track T c , is also found around a/a — 0.96. If Tj does 
in fact track T c , then this saturation may be connected 
to the filled spin up 4/ and empty 5d bands starting to 
overlap causing a mixed configuration of the J=7/2 mag- 
netic Eu 4/ 7 configuration and the J=0 non-magnetic Eu 
4/ 6 configuration. [53] Within our calculations for EuO 
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with lattice parameters smaller than the IMT values, a 
J=7/2 magnetic configuration is always generated as the 
real many-body nature of the J=0 singlet state is not 
accessible from the LDA framework. Therefore, our cal- 
culations do not include the possible disruption of the 
magnetic order from the / 6 configuration and we find T c 
continuing to rapidly rise after the gap has closed. As the 
material becomes more metallic with increasing doping, 
other exchange mechanisms such as RKKY would start 
playing an increasingly important role. The transition, 
however, will be smooth as the Fermi wave vector, kf, 
in the metallic state will be very small due to the low 
doping level. 

The mean field T c for the biaxial stress case is also 
included in Figure [2j>. This is calculated as T c = 
(2/3)S(S + l)((4jf + 8Ji) + (4j| + 2j£)). The in- 
plane (||) and out-of-plane (_L) nearest and next-nearest 
neighbor distances are calculated from the optimized FM 
tetragonal structure at a given a lattice parameter, and 
the appropriate J\ and J2 exchange parameter values are 
determined for these distances from the data presented 
in Figure [2^i. The biaxial stress shows a smaller effect on 
T c for a given in-plane lattice parameter change as com- 
pared to the isotropic stress case. This occurs because 
the average atomic distances change more slowly in the 
biaxial stress state due to the ability of the out-of-plane 
lattice spacing to compensate for changes in the in-plane 
lattice parameter. The IMT is now found to occur at 
a/a = 0.94. As in the isotropic stress calculation, we 
find T c continuing to increase as lattice parameters are 
decreased beyond the IMT. 

The compensation for the in-plane lattice parameter 
change in the biaxial stress case is shown in the inset of 
Figure [2Jl. The relationship between the in-plane, a, and 
out-of-plane, c, lattice parameters for the biaxial stress 
state is given by = = , where v is poisson's ratio. 

v is calculated to be 0.15. This is very close to the 
experimental value of 0.2 for NiO,[27 and the calculated 
value of 0.2 for GdN.gg 

Magnetism in EuO is related to three main types of ex- 
change mechanisms: local, or on-site exchange between 
the Eu 4/ and other Eu orbitals; a Kramers-Anderson 
superexchange 28J; and a number of other virtual ex- 
change possibilities across the semi-conducting gap. As 
expected we find that on-site exchange between Eu 4/ 
and 5d is significant. It will not, by itself, generate long 
range FM order, but plays a major role in the various vir- 
tual exchange mechanisms. Another conceptually impor- 
tant mechanism is the Kramers- Anderson superexchange 
caused by the overlap of Eu 4/ and O 2p orbitals. How- 
ever, since the O 2p-Eu 4/ hopping parameters are quite 
small, as compared to O 2p-3d hoping in transition metal 
oxides, this antiferomagnetic component of J 2 plays only 
a minor role in EuO. If the gap between the O 2p-Eu 
4/ states were to drastically reduce, increasing the over- 
lap in energy of the states, this mechanism would be- 
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FIG. 3: The mean field T c as a function of electron doping for 
bulk (a = c = do = 5.158A), isotropic stress (a = c — 4.92A) 
and biaxial stress (a = 4.87A, c = 5.3lA). T c of Gd doping 
from Ref. El 



come much more pronounced. Nonetheless, this can only 
be achieved by changing the chemistry of the europium 
chalcogenides; lattice parameter changes are not enough. 
Although the exact decomposition of all possible virtual 
exchange mechanisms due to the excitations across the 
semiconducting gap is outside the scope of this work, by 
using a constrained potential in LSDA, one can change 
the energy of a specific orbital and effectively control the 
hopping integrals associated with that orbital. This helps 
to separate out what types of changes in the band struc- 
ture are most significant to the ferromagnetism. 

We find that in order to most effectively increase T c , 
the hybridization between Eu 4/ and bd should be in- 
creased, as should the that of the O 2p and Eu 5d bands, 
while at the same time minimizing the O 2p to Eu 4/ 
hybridization. Since the Eu 4/ band sits between the 
O 2p and the Eu 5d in EuO it is not possible to do all 
three at the same time. However, the gap between the 
Eu 4/ and 5d bands closes much more rapidly than the 
O 2p and Eu 4/ band gap, so there is minimal lose in the 
ferromagnetic exchange due to the O 2p-Eu 4/ superex- 
change from decreasing lattice parameters. Therefore, 
we find that decreasing lattice parameters in EuO is an 
ideal way of strengthening the ferromagnetic exchange. 

Finally, by looking at the total energy changes due to 
a rigid band shift of the chemical potential for the FM, 
AFMj and AFM/j configurations, we can also study the 
combined effect of isotropic or biaxial stress and elec- 
tron doping on the mean field T c . Figure [3] shows these 
results for three situations; bulk lattice parameters, an 
isotropic stress case on the insulating side of the IMT 
(a/ao — 0.955), a biaxial stress case on the insulating 
side of the IMT (a/a = 0.945). At low doping, T c im- 
mediately increases in all cases because of the gain in 
total energy of the FM configuration relative to either 
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of the AFM configurations. As doping increases, the en- 
ergy cost to add extra electrons into both FM and AFM 
configurations becomes similar and a maximum in T c is 
reached. Our calculations suggest that the DOS of a 
5.5% biaxially compressed epitaxial film would allow a 
significant increase in T c , from 138 K to 175K. 

Also included in Figure [3] for comparison, are data 
from Gd doping experiments of Ott et al. [B] In essence, 
Gd doping is adding one electron per Gd atom to EuO 
while leaving the S=7/2 spin configuration unchanged. 
The maximum change in T c with Gd doping is about 
twice the size expected for electron doping via a rigid 
band shift of the chemical potential. This suggests that 
there are some further effects present. One possibility is 
that subtle changes occur in the 5d6s bands which effect 
the 4/ to 5d exchange. The presence of the smaller Gd 
3+ ion may also provide a positive chemical pressure on 
the EuO matrix and decrease lattice parameters, activat- 
ing all the previously discussed mechanisms. It is worth 
noting that if the impressive T c enhancement from Gd is 
related to chemical pressure, then Gd would be a less ef- 
fective means of increasing the T c of a film that is already 
biaxially compressed. 

In conclusion, we find that epitaxial influence can be 
used to increase the T c of EuO. However, we show that 
the effect of the biaxial stress generated by epitaxy is 
not as effective as hydrostatic pressure due to the ability 
of the out-of-plane lattice parameter to partially com- 
pensate for the in-plane lattice changes. The amount 
of out-of-plane compensation is minimized for the three- 
dimensional isotropic rock-salt-type structure in EuO be- 
cause of its small poisson's ratio. A material with a two- 
dimensional electronic structure would allow for the ulti- 
mate minimization of the out-of-plane compensation in- 
herent in the biaxial stress state of epitaxy, although the 
likelihood of buckling or other structural modifications 
would be enhanced. From the band structure perspective 
we find that the location of the near Ef bands suggests 
that a decrease in lattice parameters is an ideal means 
of increasing T c , as it strongly enhances the main ferro- 
magnetic exchange mechanisms. The closing of the semi- 
conductor gap is found to occur at a 5% lattice reduction 
in the isotropic stress case, or 6% reduction in the biaxial 
stress state. Although the closing of the gap will change 
the next nearest neighbor exchange mechanism, there is 
no sudden change in the magnetic properties. Finally, we 
find that electron doping will lead to an increase in T c 
for isotropically and biaxially compressed films. 
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